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Abstract 

We argue that the order of the chiral transition for iVj = 2 is a sensitive probe of the QCD 
vacuum, in particular of the mechanism of color confinement. A strategy is developed to investigate 
the order of the transition by use of finite size scaling analysis. An in-depth numerical investigation 
is performed with staggered fermions on lattices with Nt = d and Ng = 12,16, 20, 24,32 and quark 
masses amq ranging from 0.01335 to 0.307036. The specific heat and a number of susceptibilities 
are measured and compared with the expectations of an 0(4) second order and of a first order 
phase transition. A second order transition in the 0(4) and 0(2) universality classes are excluded. 
Substantial evidence emerges for a first order transition. A detailed comparison with previous 
works is performed. 
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I. INTRODUCTION 


Nf = 2 QCD can provide fundamental insight into the mechanism of conhnement. A 
schematic view of the phase diagram is shown in Fig. Q The quark masses are assumed 
to be equal for the sake of simplicity: = m; /r is the baryon chemical potential. 

Consider the plane /r = 0. As m — cx) quarks decouple and the system tends to the 
quenched limit. There the deconhning transition is well understood: the transition is an 
order-disorder hrst order phase transition, the symmetry involved is Z 3 and the Polyakov 
line {L) is an order parameter. In the presence of quarks Z 3 is explicitely broken and (L) is 
not a good order parameter. Empirically, however, it works as an order parameter at quarks 
masses down to m ~ 2.5 — 3 GeV. 

At m ~ 0 a chiral phase transition takes place at Tc ~ 170 MeV, from the low temperature 
phase where chiral symmetry is spontaneously broken to a phase in which it is restored: the 
chiral condensate ('0'0) is the corresponding order parameter. At some temperature T 4 > 
also the Ua{^) symmetry, which is broken by the anomaly, is expected to be effectively 
restored. 

It is not understood what the chiral transition has to do with the deconhning transi¬ 
tion, but empirically the Polyakov line has a rapid increase at the transition temperature, 
indicating deconhnement. 

More generally the transition line in Fig. ^ is dehned by the maxima of a number of 



FIG. 1: Schematic phase diagram of Nf = 2 QCD. 
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susceptibilities (Cy, Xm, • • •) which all coincide within errors, and which indicate a rapid 
variation of the corresponding parameters across the line. 

A renormalization group analysis plus e-expansion techniques can be made at m ~ 0, 
assuming that the relevant degrees of freedom for the chiral transition are scalar and pseu¬ 
doscalar helds 0 , 0 , Ql , or more precisely that the order parameters are the vacuum expec¬ 
tation values (v.e.v.) of the following helds 


0 : (l)ij = {qi{l + X5)qj) ii,j = . 


( 1 ) 


Under chiral and U^(l) transformations of the group 0 transforms 

as 

0 ^ ( 2 ) 


so that by the usual symmetry arguments, and neglecting irrelevant terms 




TT 


5(i(Tr{0V}) -Y^ 2 Tr{( 0 V)^}+c det 0 -h det 0 ^ 


( 3 ) 

The last term describes the anomaly: indeed it is SU{Nf) (g) SU{Nf) invariant, but not 
17^(1) invariant. 

A second order phase transition corresponds to an infrared (IR) stable hxed point. For 
Nf > 3 the effective action is of the form Eq. © and no such point exists so that the 
transition is hrst order. For Nf = 2, det(/> has mass dimension 2 so that other relevant terms 
emerge, like (det0 -|- det0^)^ and Tr{0l0}(det0 -t- det0^). If the anomaly term in Eq. 
vanishes (c ~ 0), i.e. if the r]' mass vanishes at Tc, then there is no IR stable hxed point and 
the transition is hrst order. If instead c 7 ^ 0 the symmetry is SU{2) 0 SU{2) ~ 0(4) and a 
hxed point exists which can produce a second order phase transition. 

In the hrst case the phase transition is hrst order also at m 7 ^ 0 and most likely up to 
m = 00 . 

In the second case a phase transition is only present at m = 0, which goes into a contin¬ 
uous crossover as m 7 ^ 0 : this is true also in presence of a small chemical potential /i 7 ^ 0 , so 
that a tricritical point is expected in the T-/i plane (see Fig.^) at the border of the crossover 
with the hrst order transition line which takes place in the small T, large p region 0. Pro¬ 
posals exist to detect the tricritical point in heavy ion collisions: obviously no such point 
exists if the transition at /i = 0 is hrst order. 
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The issue is in fact fundamental. If confinement is an absolute property of the QCD 
vacuum and the deconhnement transition corresponds to a change of symmetry (order - 
disorder), then a crossover is excluded and the only allowed possibility is that the transition 
is always hrst order. The argument also extends to the case of 2 + 1 flavors, of which Fig. ^ 
is a boundary j^. The question deserves a careful study. 

A few groups have investigated the problem on the lattice with staggered Q,Q,y,y, 
HHQ or Wilson fermions. The strategy used has either been to look for signs of 
discontinuity at the transition, or to study the dependence on m of the peak of different 
susceptibilities, or to study the magnetic equation of state. No clear sign of discontinuity 
has been observed, but also no conclusive agreement of scaling with 0(4) critical indexes. In 
particular the thermal exponent y* = 1/z/ (see Sect. O for the definition) determined using 
staggered fermions differs signihcantly from that of 0(4)-0(2) (no direct determination of 
the critical exponents exists for Wilson quarks). A general tendency exists however in the 
community to consider the chiral transition second order, and the line of Fig. Qa crossover. 

In the present work we have made a big numerical effort and used large lattices attempting 
to clarify the issue. Like most of the other works we use non improved Kogut-Susskind 
action, and lattices 4 x with Lg = 16, 20, 24, 32. Some scaling violations are expected 
and a more careful study with Lt = 6 and an improved action is planned in order to control 
them. 

A preliminary account has been presented at conferences Q| . The present paper contains 
more data and a full analysis. 

The paper is organized as follows: In Section |Tll we explain the strategy used to attack 
the problem. Section IIIII contains details about the simulations and the numerical results. 
Section |TVl contains the analysis of scaling. Section |V| contains a discussion and the conclu¬ 


sions. 


II. STRATEGY 

The theoretical tool to investigate the order of a phase transition is hnite size scaling Q, 
Q]. The extrapolation from hnite size Lg to the thermodynamical limit = cxo is governed 
by the critical indexes, which identify the order and the universality class of the transition. 
Approaching the transition, for a higher order or weak hrst order transition, the correla- 
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tion length of the order parameter ^ goes large compared to the lattice spacing a, so that 
the dependence of physical quantities on a/( can be neglected. More precisely, if C/kT is 
the effective action (density of free energy) 

(4) 

the dependence on a/^ disappears as Tc is approached, since ^ diverges as 




-T^OO 


r 


(5) 


where r = 1 — The variable Lg/^ can be traded with and the scaling law follows 

T ~ L-^ (ri;/-, am,Lf) . (6) 

The problem has two scales, ^ and 1/mq. The effective action depends on the order 
parameter, as dictated by the symmetry, and as r —> 0 irrelevant terms can be neglected. 
The thermodynamics is described by correlators of the order parameter, which contain in¬ 
formation on the discontinuities of the thermodynamical quantities. The most fundamental 
quantity is the specihc heat, which is always guaranteed to exhibit the correct critical be¬ 
haviour, independently of the identihcation of the correct order parameter. 

For the specihc heat the scaling law is 


Co stems from an additive renormalization Q. 

For the susceptibility x of the order parameter 0(x) 


X = d^x (O(x)O(O)) - {Oy 


the scaling law is 


(7) 


( 8 ) 




(9) 


We shall discuss the question if a subtraction is needed for x as for the specihc heat in 

sectnvni 

Analogous scaling laws can be derived for mixed susceptibilities. 

The values of the indexes characterize the transition: the values relevant to the analysis 
which follows are listed in Table HI 0(4) is the symmetry expected if the chiral transition is 
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Vt 

Vh 

V 

a 

7 

/9 

6 

0(4) 

1.336(25) 

2.487(3) 

0.748(14) 

-0.24(6) 

1.479(94) 

0.3837(69) 

4.852(24) 

0(2) 

1.496(20) 

2.485(3) 

0.668(9) 

-0.005(7) 

1.317(38) 

0.3442(20) 

4.826(12) 

MF 

3/2 

91A 

2/3 

0 

1 

1/2 

3 

V^Order 

3 

3 

1/3 

1 

1 

0 

00 


TABLE I: Critical exponents. 


second order, but it can break down to 0(2) by the lattice discretization for Kogut-Susskind 
fermions P| at non zero lattice spacing. 

The scaling law in Eq. dZD for the specihc heat is valid independent of the knowledge of 
the order parameter. The scaling law in Eq. © instead is correct only if the choice of the 
order parameter is the right one. In principle the matching between m and © can be used 
to legitimate any guess on the symmetry and on the order parameter. 

The scaling laws m and (jni) are difficult to test because they depend on two variables. 
A possible strategy is to keep one of them hxed and to study the scaling with respect to 
the other. As one can see in Tabled the index yh is the same within errors for 0(4) and 
0(2) symmetry. In order to reduce the problem to one scale, we have made a number of 
simulations at different values of Lg and arriq keeping arriqLf^ fixed and assuming yh = 2.49 
which corresponds to 0(4) or 0(2). In this way as Lg is increased, arriq 0, so that the 
inhnite volume limit corresponds to the chiral transition at arriq = 0 . 

From Eq.s m and © it follows that the maxima at constant arriqLy^ scale as 

(Ov - Oo)n,ax OC 

Xmax OC (10) 

as Ls —^ oo and their positions scale as 

= const ( 11 ) 

If 0(4) or 0(2) is the correct symmetry, the values of ajv and 7 / 1 / should be consistent 
with the corresponding values listed in Table d 

Notice that Eq.s © and © involve the long range part of the correlations, i.e. they are 
related to the infrared regime, and are not expected to be significantly affected by scaling 
violations Oiaj^. 
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If the answer to this test is positive the chiral transition is second order at arua = 0 and 
a c_ a. a™, , 0. K iaataad the aaewet is negative and the assuntption of 
about the relevant degrees of freedom is correct, the transition is hrst order at arug = 0, and 
also at anig ^ 0. 

An alternative strategy can be as follows. At hxed arUg, /3 the values of the susceptibility 
should converge at large Lg if C is analytic. One can change variable by replacing 
with the ratio 


tLV'^ 


{arngLif^y/ 

The scaling laws Eq.s dlD and m become then 


= T{amg) _ 


Cv-Co^ amgLf 

X — {T{arng)~^^^‘'y^\ arrigiy^ 


( 12 ) 

(13) 

(14) 


At large Lg the dependence on arUgLy must cancel the dependence on Lg in front of the 
scaling functions in Eq.s m and (ini). It follows that 

Cv-Co- (r(am,)-i/("2''‘)) (15) 

The peaks of {Cy — Cq) and of y should then scale as 

{Cv-Co)m.^ oc (am,)-“/(""'‘) 

Xmax OC (17) 

as arrig ^ 0. As for the position of the maxima, it scales according to 

= const. (18) 


An alternative possibility is to keep the scaling in the form of Eq.s o and 0. and require 
that the volume dependence disappears at tL^’^ hxed. This kind of scaling could work 
better if the correlation length is comparable to Lg, while aLgm.,^ 3> 1. This implies the 
scaling laws: 

Cv-Co^ {amgY^'^'^^^yc [tYJ'') 
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(19) 

( 20 ) 



Eq.s fjl7|l for the maxima stay unchanged, but the positions of the maxima scale now as 


= const 


( 21 ) 


and the width of the peaks are volume dependent. 

All that is expected to be true at sufficiently large values of aLg ■ m.^ and at sufficiently 
small values of arriq, such that we are not too far from the critical point. 

r = 1 — T/Tc is usually taken in the literature QQ as proportional to /do — /d, where 
/do is the value of /d = 2 Nc/g‘^ at the chiral {aniq = 0) transition, and all the analyses of the 
scaling law are based on that choice. 

In fact, since 

1 


T = 


the correct definition is 


Lta{/3, aniq) ’ 

T ^ a(/do,0) 


^ Tc ^ a(/d, arriq) 


( 22 ) 


(23) 


and the dependence on aniq is non trivial (see e.g. [18|). a{P,amq) is expected to be an 
analytic function in a neighborhood of the critical point and therefore for sufficiently small 
/dc — /d and m 


dd dd 

a(/d, aniq) ~ a{po, 0 ) + ^(^o, 0)(^ - /?o) + ^ 0)amg . 


(24) 


If needed higher orders in amq and (/d—/do) can be included. It then follows that at sufficiently 
small values of am„ 


with 


r = C(/do - /d + kmaniq) 


C ^ ^(/3„.0) 


(26) 


k-m, — 


df3 
1 5 In a 
C damn 


iPoy 0 ) . 


(26) 


In the quenched case this reduces to r cx /do — /d as usual: in the presence of dynamical 
quarks 7 ^ 0 |l^ . 

The scaling law for the position of the peaks becomes then 


/do — /dc + kmCimq = const ■ 


(27) 
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Eq.s (fTH|l and (fTH|l should be valid if S> ^/a, Lg 3> l/(am 7 r) (see Table |TT)). If the 
alternative possibility is considered, i.e. requesting that the free energy stays hnite when 
Lg —> cx) at hxed the position scales instead as 

(^0 — f^c + kmarriq = const ■ . (28) 


Analogous formulae can be written including quadratic terms of the expansion Eq. (1241) (see 
Section HVl below). 

An alternative technique is to investigate the order of the transition by looking for dis¬ 
continuities: if the transition is hrst order at m = 0 , it is expected to be so also at m 7 ^ 0 . 
If the transition is weak hrst order, at small volumes compared to some critical volume it 
will behave as if the free energy were regular, so that Eq.s (USD, (uni) and (HZD , or (m and 
(EOl), are expected to be valid, with the critical indexes appropriate to hrst order. At larger 
volumes, however, the peak of the specihc heat as well as the peaks of the other susceptibil¬ 
ities should increase proportionally to the volume, as a consequence of the discontinuity in 
the hrst deriva^ves of the free energy. At the same time a bistability should appear in the 


time histories [S|. Such an analysis has been in particular performed in Ref. |l lL h2||: some 


sign of growth with the volume has been observed, but no signihcant bistability; we will 
comment on this result below. Of course if a discontinuity is observed one can conclude that 
the transition is hrst order. If not one cannot exclude that it could be observed at larger 
volumes. 


Finally one can investigate the so-called magnetic equation of state 1^, i.e. the scaling 


behaviour of the chiral order parameter itself, (V’'*/’), versus the reduced temperature. The 
scaling law is 


firm 


(29) 


and again it can provide information on the critical indexes. 


III. NUMERICAL SIMULATIONS 


A. Algorithm 


Monte Carlo simulations were performed using the standard staggered action. The Hybrid 


R algorithm 


17[ was used for the conhguration updating. Since it is a non-exact algorithm. 
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its systematic errors must be kept under control. The finite integration step used in the 
molecular dynamics evolution introduces a systematic error on the mean values of observ¬ 
ables proportional to the integration step squared Ar^. Great care was taken to ensure that 
this systematic shift were much smaller than the statistical error in each Monte Carlo run. 
Typical values of the integration step At vary with the mass of the quarks in units of the 
lattice spacing as At = airtq/A. When the use of that value for the integration step was 
too proibitive, i.e. at the smallest quarks masses used in this work, the integration step 
was in any case taken below amq/2. The stopping condition used for the conjugate gradient 
inversion was fixed requiring that the residue were smaller than 10“®. The length of the 
molecular dynamics trajectories was fixed to 1 for all of our simulations. 

B. Run parameters 

All of our numerical simulations were performed using a lattice temporal extent of Lt = 4. 
To begin with we run two sets of Monte Carlo simulations fixing for each the value of arngL^^ 
as explained in Section |nj The two sets, called in the following Runl and Run2, have 
arngL^A = 74.7 and arngL^^ = 149.4 respectively. The spatial lattice sizes Lg used for each 
of the two sets are Lg = 12,16, 20, 32. 

Additional simulations at Lg = 24 and arrig = 0.04444 and at Lg = 16 and arug = 0.01335 
were added. The second one was chosen by purpose at the same mass of the Lg = 32 of 
Runl. 

A summary of the bare quark masses and Lg used is reported in Table m The total 
number of MC trajectories collected is also reported together with the quantity aLg ■ rriT^ at 
the pseudocritical value of the coupling {3c {am.^ was taken from the parametrization given 
by the MILC collaboration in Ref. Q), S.ce fo. ou. the epat.a, extent . ntnch 
larger than the pion correlation length, no large infrared cut-off effects are expected, except 
possibly for the run at Lg = 16 and arrig = 0.01335 (see Table HH). 

For each value of arrig and Lg and for each run, MC simulations were performed at 
different {3 values in order to inspect and to have under control the whole interesting critical 
region. See Appendix El for the whole listing of our run parameters. 
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C. Data Reweighting 


The collected raw data were analyzed using standard statistical procedures (see e.g. Iia). 
For the history of each observable, thermalizations were taken self consistently to be hve 
times the integrated autocorrelation time estimated from thermalized trajectories. 

The data collected were analyzed using the multi-histogram reweighting technique com¬ 
bining data taken at different /d’s together. This method allows to extract mean values of 
observables and their susceptibilities at intermediate (3 values over the whole range explored 
with numerical simulations. Using the reweighted data, it is possible to locate accurately 
the position at which the susceptibilities attain their maximum and their value at the maxi¬ 
mum. Sometimes in previous studies a single point with high statistics at about the critical 
point was used. Using data from simulations done at several (3 values covering the whole 
critical region, eliminates the risk of a wrong extrapolation from a single (3 too distant from 
the critical point^. Moreover the method allows a better sampling of the probability distri¬ 
bution, due to the fact that different simulations are combined together, thus increasing the 
precision and conhdence of the measurement. 

The errors of observable quantities were estimated using the bootstrap method. In prac¬ 
tice, this means repeating the whole multi-histogram reweighting procedure a number of 
times starting from random data samples distributed as the measured empirical distribu¬ 
tions. 


TABLE II: Run parameters for the numerical simulations. 



Runl 

Run2 

Other 

Ls 

12 

16 

20 

32 

12 

16 

20 

32 

16 

24 

aniq 

0.153518 

0.075 

0.04303 

0.01335 

0.307036 

0.15 

0.08606 

0.0267 

0.01335 

0.04444 

# Traj. 

22500 

87700 

14520 

14500 

25000 

131390 

16100 

15100 

10000 

10000 

aLg ■ rriT^ 

11.9 

11.0 

10.0 

8.9 

11.3 

15.8 

14.8 

12.4 

4.5 

12.2 


^ Remember that for single histogram reweighting the statistics needed in order to extrapolate measured 
quantities at a value of /3 = do + different from that used in the actual simulation grow exponentially 
with Ad. 
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D. Observables 


For each generated configuration of our MC simulations we measured the average spatial 
and temporal plaquettes (Pq-, Pr), the chiral condensate {'iJj'iIj), the energy density 
and the following lattice susceptibilities (the notation is the same as in ll|): 


disc _ 
X.m 




1 r 

V - 


((TrP-1)") - (TrP-^)2 


conn 


Ay ^ \ x,y 

^,y 


Xe,« = r[{FiFj) - (P.) (Fj)], i,j = a,T 

\'e,i = yl{PM’Doi’)) - {Pi} (i’Datpy, ! = ff, r 
Xt.i = Vl^Ptiinpy - (Pi) (iiip)], i = a,T 
Xtj = V[({inl>){iiD(ii>)) - ({inp}) (ipDoi>] 


( 30 ) 

( 31 ) 

( 32 ) 

( 33 ) 

( 34 ) 
( 36 ) 
( 36 ) 

where V = L^Lt is the yolume; Dq is the temporal component of the Dirac operator D. 

The connected component of the chiral susceptibility been measured for all 

of our simulations but only for a fraction of them. The method used to extract Xm"'^ is the 
yolume source method without gauge fixing as described in Ref.jll|. The disconnected com¬ 
ponent giyes the dominant contribution for large yolumes and small masses. The connected 
part is instead releyant at small yolumes and relatiyely large masses. For most of our lattices 
we haye determined the connected part only around the peak, and we haye considered it 
as a constant with respect to /?. We estimate that this is a good approximation within our 
errors. A representatiye example^ is shown in Fig. |21 (taken at Lg = 24, arug = 0.04444). 
The typical contribution to Xm of is less than about 15% of Xm^'^ fhe peak yalue 
and is a slowly yarying function of (3. 

A comprehensiye list of measured yalues of (Po-), (Pr), (V’V’)) (V’-^oV’) and the lattice 
susceptibilities for our MC simulation can be found in Appendix El 


For this lattice y™”" was measured for all points. 
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1 M 
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5.3 5.31 5.32 5.33 5.34 5.35 


FIG. 2; Comparison between the connected and disconnectd component of Xm- The former is 
typically a small fraction of the connected component at the peak position. 

IV. SCALING ANALYSIS 


The basic thermodynamic susceptibilities, i.e. the specific heat Cv, the chiral suscepti¬ 
bility Xm and the mixed susceptibility Xt 

1 


Cv = 


Xm 


Xt = 


VT^ a(l/T)2 
T 


V drri^ 
T 


InZ 




In Z (37) 

(38) 

V o{l/T)omq 

can be expressed as sums of the lattice susceptibilities O-dSSl) multiplied by regular func¬ 
tions of [3 and arug. Specifically the specific heat Cy is a function of Xe,iji Xe,i and Xej] 
Xt is a function of Xt,i and Xtj- The contribution of other susceptibilities entering in the 
expression of Cy and Xt involving the quark mass are expected to be negligible in the chiral 
limit j^. As for Xm, d/d{mg) in Ea. dH^ is intended at constant temperature. Since tem¬ 
perature depends not only on [3 but also on arrig, the physical Xm is a combination of x'^^i 
Xty and Xeyr witli computable coefficients. 

By Cy in the following we mean Xe,aa] the analysis with Xe,i and Xej is similar and 
compatible with respect to scaling. 

We are interested in studying the singular behavior of Gy, Xm and Xt as the critical 
surface is approached which is given by the most singular divergent quantity among the 
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lattice susceptibilities corresponding to a given termodynamical susceptibility. 


A. Pseudocritical coupling 


One of the observables analyzed in the literature to understand the order of the transition 
has been the position of the peaks of thermodynamic susceptibilities as a function of arug. 
The position of all these peaks happen to coincide at given arUg and Lg, thus defining a unique 
(pseudo)critical coupling Pc^aniq). Previous works in the literature assume r oc /do — /9, a 
choice usually based on a strict analogy between QCD and the 0(4) statistical model. In 
fact the correct thermodynamical reduced temperature is given by Eq. (12dj) . In principle 
the dependence of a(/3, arUg) on arUg could be measured by use of independent quantities 
(see e.g. ^^). We will try a fit of the position of (3c by a form like Eq. (EZI) or (j2Hl), which 
is expected to be valid at sufficiently small values of arUg. To extend the range of validity 
of the approximation the quadratic terms proportional to am^, amg{(3o — (3) and {(3q — (3)“^ 
may be added: 


r cx (/?o - /?) + kmarrig + k^ 2 {amgf + kmf}amg{(3o - (3). (40) 

A term kjj2{(3Q — /3)^ turns out to be negligible. In fact one can write the lattice spacing a 
as: 

^1 

where the deviation from asymptotic scaling are represented by the fact that the term 
A(/d, arUg) is [3 dependent. This term is slowly varying with (3 and can be well described 
by a linear function of (3 with coefficients depending on arUg in the relevant range of /d’s. 
The fc /32 coefficient is given by lna/df3'^\p=i3g^amg=o can be computed and can be neglected 
within errors. The other unknown parameters can be fitted to the data. 

Fig. H shows the critical line (3c{amg). Our determinations are reported together with a 
collection of world data^. A good agreement among different determinations can be appre¬ 
ciated. 


^ Data of the JLQCD collaboration are taken from Ref. lul . We thank E. Laermann and C. DeTar for 
providing us with the data of the Bielefeld group and of the MILC collaboration respectively. 
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FIG. 3: (Pseudo)critical couplings determined in this work are shown together with other de¬ 
terminations form the literature. Logarithmic scale is used for the horizontal scale in the main 
figure, while a linear scale is used for the smaller inset figure. The dotted line is the best fit 
(x^/d.o.f = 0.89, d.o.f. = 15) curve including arrig and amg(/3o — Pc) terms for masses arug < 0.4. 

The expected variation of r as a function of anig or Lg is given by Eq. (nm or (EH): 

r = or (42) 

r = KLlI- . (43) 

Notice that for a first order transition the exponent of arUg is 1 and the term on the right 
hand side of Eq. (H2|) can be reassorbed in km so that such term can be discriminated only 
for a second order scaling behavior. 

The unknown coefficients km, km?, kmp and kj. of the expansion of r around the critical 
point are determined by use of a best fit procedure. We start from the form Eq. (liHj) . In 
agreement with previous works we find that the position of the peaks does not depend on 


the lattice size 


I. llTl. i.e. that fc). = 0 within errors implying that no information can be 
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obtained about the order of the transition (first order or second order 0(4) or 0(2)). The 
quality of the £t assuming /c). = 0 is shown in Fig. Eland the resulting coefficients are listed 
in the first line of Table nm These coefficients are obtained by a best £t up to a maximum 
value of arrig, {amg)max, which is then extrapolated to zero. They are stable and consistent 
with a linear fit at low values of arug (< 0.0267) as shown in the second line of Table HTTl 

The scaling of Eq. m or assumes that aLgrriT, ^ 1 but the correlation length ^ can 
be comparable with Lg, which is certainly true sufficiently close to the critical point in case 
of a second order or weak first order chiral transition. 

We have then analyzed the scaling of the form Eq. m or If the transition is first 
order the analysis coincides with the analysis done for the scaling Eq. (P|) and km is in fact 
h — k 

If the transition is 0(4) a similar analysis can be performed (similar results hold for 0(2) 
or mean field). The y^/d.o./ is acceptable, the result for the coefficients is shown in the 
third line of Table Uni kr is consistent with zero, and the result is therefore compatible 
with that of the first line. However the £t becomes unstable if we try to extrapolate to low 
masses keeping only the linear term of Eq. (j4()jl (line 4 of Table nTT|l . 

The critical coupling Pq and the coefficient km are stable both for first order and 0(4) 
behavior and can thus be conhdently estimated. Our final estimates, obtained by a weighted 
average of linear and quadratic hts, are Po = 5.2484(5) and km = 1.82(8) for a first order 
transition; Pq = 5.2435(25) and km = 1.13(19) for 0(4). Other terms cannot be reliably 
estimated with present data. In particular we are not able to discriminate the contribution 
of quadratic terms from that coming from kr and consequently it is not possible to establish 
the order of the transition by looking at the (pseudo)critical couplings alone. 

TABLE III: Best fit parameters for the scaling of the pseudocritical coupling. Different kind of fits 
are explained in the text, k'^ is set to zero as explained in the text. 


Po 

km 

kr 


km/3 

5.2484(4) 

1.84(7) 

= 0 

-0.3(2.4) 

-4.3(2.7) 

5.2481(13) 

1.75(13) 

= 0 

= 0 

= 0 

5.2430(40) 

1.20(60) 

-0.1(1) 

-7(6) 

6(6) 

5.2437(21) 

1.12(16) 

-0.134(46) 

= 0 

= 0 
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FIG. 4: Different behavior of predicted (pseudo)critical couplings for a first order transition and a 
second order 0(4). Present data do not permit a clear discrimination using only jdc- The lu band 
is displayed. Continuous lines correspond to a first order transition while the dotted lines are the 
prediction for the 0(4) symmetry. 

The possibility to discriminate between the hrst order and 0(4) behavior from the mea¬ 
surement of f3c is in practice very faint. Fig. (j^ shows the different predictions for the 
(pseudo)critical coupling based on available data. If we discard the high values of the masses 
a possible difference would only be visible at very small bare quark masses. Using estimates 
of kr shown in Table llfll a quark mass of about 0.003 should be used, which requires a big 
numerical effort. 

We would like to remark that the explicit dependence of r on arug is in any case necessary 
to fit present data. Fitting to a function of the form 

/?c = ^0 + Cr{amqy*^y^ (44) 

with 0(4) values for the exponents, gives a x^/d-o.f ~ 21 in the interval up to {amq)max = 
0.075, which decreases as {amq)max decreases and is ~ 2 at {amq)max = 0.02. 
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As a final remark, the dependence on (3 and amq of the lattice spacing a Eq. |2S1 can be 
measured from other observables (see e.g. In particular our estimate for km for a hrst 

order transition (line 1 of Table llll|l is compatible with those of Ref. km ~ 1-95 (affected 
by errors of order 20%). We notice that = 0 implies that r = 0 on the critical line, or, 
by Eq. (El that the critical temperature is independent of amq near the chiral point. 

B. Scaling at fixed amqL'^’^ 

As explained in details in Section IHl we have adopted a novel strategy in order to simplify 
the two scales problem. We assume 0(4) — or 0(2) — critical behaviour and we use this 
assumption to £x a dependence between amq and Lg in our Runl and Run2 so as to £x the 
second scaling variable in Eq. and reduce the problem to a one scale problem: in this 
case the only assumption is 0(4) itself. This allows us to test whether 0(4) is consistent or 
not with data without any further approximation. 

We hxed amqL^^ = const with yh = 2.49 that is the value expected for 0(4) and 0(2) 
critical behavior, with const = 74.7 for our Runl and const = 149.4 for Run2. The following 
scaling formulas should hold (see Eq.s m and Q): 


Cy(r , Lg) — Cq 

= L':/''<s>c(.TLi/n 

(45) 



(46) 


and in particular the peaks of the specihc heat and Xm should scale as Eq.s (HDD. 

The subtraction of the non critical part Cq for the specihc heat is needed. In principle 
it can be obtained as a parameter from the £t to the maxima of the specihc heat. However 
since we also have data at /3’s diherent from the pseudocritical coupling, we were able to 
perform a direct measurement of this quantity. Appendix ^reports the details of the study 
of the background for Cy- Our hnal estimate for the background is C'o(/9) = 0.400(43) — 
0.0663(83)/3. Note that the /3 dependence is very weak and assuming a constant value for 
the background Cq does not modify the following analysis. No dependence of Cq on amq is 
observed. 

For the susceptibility of the chiral condensate Xm, for the moment we do not operate any 
bakground subtraction. 

The measured peak values for the subtracted specihc heat Cy — Cq and chiral condensate 
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FIG. 5: Specific heat (top) and Xm (bottom) peak value for Runl (left) and for Run2 (right), 
divided by the appropriate powers of Lg (Ea.s H^lHHl) to give a constant. Both the 0(4) and 0(2) 
critical behaviors are displayed. Notice that for the case of Xm the ratio 7 / 1 / have almost the same 
numerical value so that the two curves are almost indistinguishable. 

susceptibility Xm for Runl and Run2 are shown in Fig. El They are evaluated on the curve 
obtained by reweighting. The figure shows the peak values of susceptibilities rescaled by the 
appropriate power of the spatial lattice size Eq. (HUD- If the scaling laws m would hold, 
the displayed quantity should be a constant. Visibly this is not the case. 

The 0(4) and 0(2) critical behavior is clearly in contradiction with the lattice observation. 
In particular 0(4) and 0(2) scaling predicts no singular behavior in the Lg ^ 00 limit for 
the specific heat as the critical exponent a is negative. This means that as Lg is increased, 
the singular part of Cy should decrease with volume, i.e. that the specihc heat should not 
grow which is in clear contrast with the data. Also for the chiral condensate susceptibility 
Xm the predicted exponents fail to reproduce lattice data. In either case the x^/d-o.f. of the 
£t with a constant function excludes the behavior of scaling law m- The full scaling laws 
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FIG. 6: Scaling of the specific heat (top) and Xm (bottom) for Runl (left) and for Run2 (right), 
see Eq.s 0 and ©• The curves are obtained by reweighting. 

Eq.s (0) and (©) were also studied (see Fig.(jn))). The horizontal scale was obtained by htting 
the pseudocritical temperature as described above. As one can expect from the previous 
discussion, data don’t scale according to the predicted laws. 

Similar hgures are obtained assuming 0(2) symmetry. 


C. Scaling at Lg 


oo 


A further scaling test can be done supposing that the lattice size is much larger that 
all other relevant physical lengths. In such a case one expects that the system show the 
same behavior of an inhnite system. The scaling laws expected in this case are those of 
Eq.s fTHjl and (fTT)|) . These equations predict no dependence on the parameter L.. This is 
the same assumption used in previous scaling analyses of the chiral transition (iCi IllL Il2||. 
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FIG. 7: Specific heat (top) and Xm (bottom) peak scaling for 0(4) (left) and first order (right). 

The alternative possibility, illustrated in section ITTl is to keep (i.e. ^/{aLg)) fixed 

thus remaining with the scaling equations Eq.s dH and (EDI). Physically this means that 
the correlation length ^ is not small compared to Lg, which is certainly true in the vicinity 
of the critical point. 

These scaling laws are only expected to hold for large values of and small masses. 

The difference between the two alternatives is only visible by considering the width of 
the susceptibilities peaks, the heights having the same behavior (see Eq. sdlZD). 

We have thus hrst performed the analysis of the scaling of the maxima of the specihc 
heat Cv and Xm- 

We have tested the different second order critical behaviors compatible with the scenario 
of Ref. Pi, namely 0(4), 0(2), Mean Field and hrst order. The peak value of the specihc 
heat and of Xm divided by the appropriate power of the quark mass should be a constant. 
These ratios are shown in Fig. m both for 0(4) and hrst order. The hgure shows also the 
conhdence region from a ht with a constant value together with the corresponding y^/d.o.f. 
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FIG. 8: Comparison of specific heat (left) and Xm (right) scaling for 0(4). Eq^l EOl (top) and 
Eql3ini (bottom). 


From the values of the x^/d.o.f. it is easily seen that the second order critical behavior is 
not compatible with data. It must be noticed that, although the validity region of the scaling 
law is not known a priori, so that the upper mass limit for the hts is somewhat arbitrary, if 
we further restrict the mass region, the x'^/d.o.f. for the Cv fits tend to increase. We stress 
that also .a ptevious studies fl Q ^ the values fouud fot the suseeptih.hty peahs wete 
not compatible with the critical indexes of 0(4), 0(2) and Mean Field. On the contrary the 
first order behavior looks compatible with data (even if the y^/d.o./. is ~ 2 ) for the specific 
heat. We will discuss the scaling of Xm in Sect llVDl 

The scaling of susceptibilities at all (3 values can also be investigated. The situation is 
depicted in Fig. |Hlfor 0(4): all the curves should coincide within errors if there were 0(4) 
scaling. Similar features are observed for 0(2) and Mean Field. No scaling is observed: this 
is clearly the case for the maxima of the susceptibilities as discussed above, but it is also 
true for the width of these curves. The analogous figure using first order indexes is shown 
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FIG. 9: Comparison of specific heat (top) and Xm (bottom) scaling for first order. EalT!llandOni 

in Fig. El Scaling is observed for the specific heat in the form Eq.s (11911 and (EOl) and not in 
the form m and (Cni), which does not describe the widths of the peaks. Again for Xm we 
postpone the discussion to Sect llVDl 



FIG. 10: Comparison between Xe,aa (left) and Xm (right) for aniq = 0.01335 and Ls = 16,32. 
Runs with higher values of the masses {arriq > 0.075) do not obey the scaling laws. The 
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one with Lg = 16 and aniq = 0.01335, which should coincide with the Lg = 32 at the same 
bare mass, in case of scaling Eq.s m and (fT^ is instead different (see Fig. cni). For scaling 
Fq.s (I10|l and (EDI), the maximum should be the same, and the widths should differ by a 
factor of 8, and this is not the case. We have carefully checked the stability of the curves 
obtained by reweighting against variations of the statistics, e.g. by discarding single data 
points from the analysis. This lack of scaling can be interpreted as due to the small value 
(4.5) of the parameter aLgUiT,, invalidating the limit bringing to Eq.s (US}, (fTK} . (ITH} and 
( 12011 . A similar effect could be responsible for the observed increase of the peak with the 

n 

volume observed in previous studies |lll |. 

A more careful study of this effect will be done in the future. 


D. Magnetic equation of state 


As a further test of scaling we check the equation of state Eq. dSH}. No scaling whatsoever 
is observed, neither 0(4)-0(2) nor Erst order, if the raw measured data are introduced in 
Eq.(j2n} (see Fig ITT} . Indeed as is different from 0 at large /3 a subtraction is needed: 
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FIG. 11: Magnetic equation of state using non subtracted data for 0(4) (left) and first order 
(right). 


the critical part of the chiral condensate has to be zero far above the critical region. A 
tentative way to understand the non critical background can be to assume it equal to the 
value {'0'0)oo of at /d = cxo. This can be computed analytically and numerically on the 
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flat configuration U^i, = 1. The result is 


{iitp) oo{amq) = 


E 


arrin 


(47) 


+ Si.i,2,3(sin 13 -)" + (sin 
which is almost a linear function in the mass range am^ < 0 . 1 . 

We then plot the subtracted value {'4>'4>)s = ii’i’) — ('0'0)oo rescaled as {^ip)s/{amqY^^ 
versus T{amq)~^^^'''^*^'> to test the scaling Eq. (El- Fig.ini shows the result for 0(4) (similar 
hgures being obtained with 0(2) and mean held). Visib 


analogous investigation has been performed in Ref. 


y 0(4) scaling is not obeyed. An 
12 j |. without the subtraction of the 
j3 = 00 term: also in that case results were in disagreement with 0(4) scaling. 

An alternative procedure is to subtract for each {aniq, Lg) the value at the largest mea¬ 
sured f3- The result is consistent with Fig. [H A third possibility is to impose scaling at 
r = 0 and look if it is obeyed at r 7 ^ 0. The result is shown in Fig. Eland again there is no 
scaling. 

If the analysis of is correct, our results, which exclude a second order critical behavior. 
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FIG. 12: Equation of state for 0(4), obtained by subtraction of ('i/’V’) at /3 = 00 . 
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FIG. 13: Equation of state for 0(4), obtained by imposing scaling at r = 0. 

should then imply a hrst order chiral transition. 

As a test of this possibility we have repeated the scaling analysis of the equation of state 
using hrst order critical indexes: the three procedures described above are consistent with 
each other and the result is shown in Fig. EJ A good scaling is observed. 

We can investigate the consequence of the subtraction needed to isolate the critical part 
of (V’V’) on the scaling of Xm- If (V’V’) — F{t /by differentiating with 

respect to aniq at hxed temperature we get 

Xm - - ^am\/^~^{T/am\/''y^)F' (48) 

oaniq 0 ^ vyh 

Keeping in mind 1/5 = {d- yh)/yh, 7 = {‘^Vh - d)/yt and v = l/yt we hnd, at r = 0 

The quantity which scales is not but a term must be 

added to it to get scaling. The subtraction of (V’)/’), due to the explicit breaking of chiral 
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symmetry at arriq ^ 0, implies a subtraction for Xm (which in the mass range of interest is 
almost constant). 

This suggests to repeat the test of scaling for Xm (Fig-s 151(1171 and |H1) by introducing a 
subtraction by a constant to be determined. The content of Fig.^Jinstead stays unchanged 
since the curves refer to the same value of the mass. 

The best £t to determine the background is done by requesting scaling for the peaks of 
Xm- the result is shown in Fig.s ED and |7D For 0(4) —and 0(2)— no scaling is obtained. 
A reasonable scaling results for hrst order. The corresponding modihed scaling at all /3’s is 
shown in Fig.s ED and n 

We notice that, due to the dependence of the temperature T on ara^ and not only on /3 the 
dehnition of ('0-0) = and Xm = must be revised, resulting in a combination 

of several terms analogous to what happens for the specihc heat. 

Analogously to what has been done for the specihc heat, we have not rehned consequently 
our analysis. Our investigation is exploratory: our Lt is not big enough, our action is not 
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FIG. 14: Equation of state for first order. 
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FIG.IHI: Scaling of the maxima of Xm — Xo at fixed amqL^^ for Runl (left) and Run2 (right). Two 


curves are shown corresponding to 0(4) (circles) and 0(2) (squares) critical behavior. The value 


of Xo is obtained by a best fit procedure. 
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FIG. 13: Scaling of the maxima of Xm — Xo at small aniq for 0(4) (left) and first order (right).The 
value of Xo is obtained by a best fit procedure. The x^/d.o.f of the best fit is also shown. 


improved. We will extensively use the correct definitions in the planned improved version 
of this investigation. 


E. Metastabilities 

A hrst order chiral transition implies hrst order also at m 7^ 0 and this should be visible 
in time histories. Metastable states should be present and should be visible as double peak 
structures in the histograms of distributions of the value of observables for large enough 
volumes. Although the results for scaling presented in the previous sections indicate that 



















FIG-ini: Scaling of subtracted chiral condensate Xm — Xo for Runl (left) and Run2 (right). The 
background value xo is obtained from a best fit of the peak value. 
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FIG. IHI: Scaling of subtracted chiral condensate Xm — Xo at small arug for 0(4) (left) and first 
order (right). The background value xo is obtained from a best fit of the peak value. 


our volumes are not large enough, we have analyzed the probability distribution function of 
a number of observables, in particular of the spatial plaquettes. The results are shown in 
Fig. nil As in previous works no convincing metastability appears. 

To estimate the probability distribution function (PDF) from a given hnite data set a 
non-zero width of the integration region should be chosen. We assume for the width a 
value given by one hundredth of the difference between the largest and smallest entry in our 
dataset. Errors are estimated using the bootstrap technique. 


29 



























































































O 50 lOO 150 200 250 

U ' I 


50 lOO 150 200 250 


FIG. 15: Probability distribution function of average plaquette. Chiral condensate shows a similar 
behavior. 

V. DISCUSSION AND CONCLUSIONS 


The study of the nature of the chiral transition of Nf = 2 QCD is of fundamental 
importance: a second order phase transition would mean crossover at m 7 ^ 0 and also at 
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finite chemical potential, implying the presence of a tricritical point in the T — /i plane 
detectable by heavy ion experiments; a hrst order phase transition would drastically change 
this scenario. Moreover this is relevant to understanding what conhned and deconhned really 
means in Nature, i.e. if they are two different phases of matter corresponding to different 
realizations of some symmetry with an order parameter, or are connected by a crossover. 

Previous studies on the subject did not come to a dehnite conclusion, mainly because of 
the huge computer power required. We have approached the problem by dedicating a big 
amount of computer power and by proposing a novel strategy for the scaling study around 
the chiral critical point. In particular we have developed a scaling study which assumes the 
critical indexes of the expected second order universality class (0(4) or 0(2) according to 


the analysis of Ref. j^) to reduce to a one scale problem without any further assumption. 
In this way we have been able to dehnitively rule out the possibility of 0(4) or 0(2) critical 
indexes. 

We have introduced the mass dependence of the reduced temperature r lEq. 1^ . which 
was neglected in previous works. 

We have then analyzed our data to test the scaling as done in previous papers, assuming 
to be already in the thermodynamical limit (Eq. fTHl and . The arrig dependence of the 
pseudocritical (3 can not discriminate between 0(4), 0(2) and hrst order. However the 
behavior of the peak of both the specihc heat and of the chiral susceptibility dehnitely 
excludes 0(4) and 0(2), but are qualitatively consistent with hrst order (Fig. [3 and ED). 

As for the shape of the critical peaks again 0(4) and 0(2) are dehnitely ruled out. The 
dependence of the width on Lg shows that the thermodynamic limit is not reached. Instead 
a scaling at Lg/^ hxed agrees with hrst order, and again dehnitely excludes 0(4) and 0(2). 

The magnetic equation of state is also nicely compatible with hrst order. 

No clear signal of metastability, except possibly some hint, is observed. 

In conclusion in = 2 QCD the chiral phase transition is not in the universality class 
of 0(4). Data strongly indicate a hrst order phase transition. Further study is needed to 
put this statement on a hrmer basis. First we are planning a similar analysis with improved 
action and L* = 6. A consistency check will also be the study of the r]' mass at the deconhning 


transition, in accordance with the analysis of j2|, and we are also working at it. 
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APPENDIX A: MONTE CARLO PARAMETERS AND RAW DATA 


In this appendix we report the details of our numerical results. Time histories of the 
observables are also available at request. 

In the early stages of the computations some of the susceptibilities were not measured. 
Blank entries in the following tables refer to missing data. 

In some previous report the data for the plaquette susceptibilities were defined with an 
extra factor of 1/4. Below we have eliminated it in agreement with dehnition Eq. to (jdbj) . 

The connected component has been measured only at the (3 value nearest to the 

critical coupling (except for the run with Lg = 24, aniq = 0.04444). 

TABLE IV: Mean value of plaquettes, chiral condensate and energy density from MC simulation. 
The number of thermalized trajectories is also shown together with the extimated autocorrelation 
time (for Pg and {ipip)) 


b 

# dVaj. 

Pa 

Pr 

ill'll)) 

(V’DoP) 

'^int 

Tint (('0V')) 


II 

Lg=VZ 


amg=[).\ 

535TE 

5.3 

1300 

0.47880(16 


0.47893(16 


1.07450(60) 

— 

7.7(2.4) 

1.47(20 

) 

5.325 

1300 

0.48578(17 


0.48601(17 


1.04638(59) 

— 

7.3(2.2) 

1.71(26 

) 

5.35 

1300 

0.49337(21 


0.49363(20 


1.01300(79) 

— 

10(3.7) 

3.90(87 

) 

5.375 

2930 

0.50082(12 


0.50130(13 


0.97909(55) 

0.72133(24) 

9.2(2.0) 

3.71(53 

) 

5.375 

1250 

0.50147(15 


0.50177(15 


0.97683(56) 

— 

6.0(1.7) 

1.91(31 

) 

5.3875 

2860 

0.50637(16 


0.50689(16 


0.95072(81) 

0.72657(24) 

12(3.2) 

6.4(1.2) 

5.3875 

1200 

0.50653(30 


0.50711(32 


0.9497(20) 

— 

18(9.3) 

15(6-9) 

5.4 

2920 

0.51129(16 


0.51196(18 


0.9233(11) 

0.73140(25) 

9.7(2.2) 

8.3(1.8) 

5.4 

1200 

0.51041(27 


0.51105(34 


0.9275(20) 

— 

13(5.9) 

16(7.7) 

5.40625 

2890 

0.51494(30 


0.51576(36 


0.8990(26) 

0.73594(34) 

32(14) 

33(14) 


5.40625 

600 

0.51464(37 


0.51582(34 


0.8986(31) 

— 

14(8.9) 

16(10) 


5.4125 

2890 

0.51762(50 


0.51869(57 


0.8796(43) 

0.74062(55) 

72(46) 

77(51) 


5.4125 

1200 

0.51955(36 


0.52100(49 


0.8635(37) 

— 

22(12) 

36(25) 


5.41875 

2900 

0.52200(36 


0.52357(40 


0.8470(28) 

0.74866(41) 

46(23) 

42(20) 


5.41875 

840 

0.52247(49 


0.52409(53 


0.8468(37) 

— 

20(13) 

28(17) 


5.425 

2970 

0.52508(18 


0.52692(23 


0.8255(17) 

0.75432(29) 

13(3.5) 

19(6.3) 

5.425 

1200 

0.52576(35 


0.52774(37 


0.8204(29) 

— 

19(9.8) 

20(11) 

5.43125 

2910 

0.52802(19 


0.52990(21 


0.8063(15) 

0.75839(24) 

17(5.4) 

16(5.2) 

5.43125 

540 

0.52952(58 


0.53151(70 


0.7917(61) 

— 

40(44) 

66(90) 

5.4375 

1300 

0.53152(17 


0.53381(16 


0.7782(11) 

— 

8.2(2.6) 

8-6(2.8) 

5.45 

1300 

0.53475(20 


0.53730(22 


0.7598(16) 

— 

10(3.9) 

12(5.1) 


32 







































TABLE IV: (continued) 


p 

# Traj. 

5.475 

1300 

5.5 

1300 

5.525 

1350 

5.55 

1350 

5.575 

1350 

5.6 

1350 


Pa 


0.54153 


0.54614 


0.55052 


0.554344 89 


0.558004 86 


0.561446 72 


Pr 


0.54411 13 


0.54917 11 


0.55340 11 


0.557310 76 


0.560962 80 


0.564575 86 




0.72134 


0.69862 


0.67937 


0.66412 


0.65103 40 


0.63826 36 


{'lpDoli>) {Tint [Pa 



4.19 95 


3.13 62 


3.03 60 


2.32 41 


Tint ((V’V')) 


5.7 1.5 


3.42 


3.41 


2.22 


2.14 35 


2.20 36 


5.3 

670 

5.33 

4840 

5.3325 

3600 

5.335 

4000 

5.3375 

4000 

5.34 

4000 

5.345 

4000 

5.3475 

4000 

5.35 

7130 

5.3525 

14800 

5.355 

4000 

5.3575 

3800 

5.36 

13200 

5.365 

4000 

5.37 

4000 

5.38 

3100 

5.4 

1200 

5.45 

1300 

5.5 

1300 

5.6 

775 


0.48822 


0.499282 83 


0.50022 


0.50164 


0.50321 


0.50425 


0.50778 


0.51042 22 


0.51126 


0.51397 16 


0.51614 


0.51776 19 


0.51941 


0.52181 


0.52429 


0.527678 78 


0.532769 91 


0.542526 64 


0.550574 49 


0.564279 58 


0.48859 20 


0.499991 91 


0.50106 13 


0.50252 14 


0.50418 14 


0.50530 24 


0.50905 23 


0.51203 26 


0.51291 21 


0.51585 19 


0.51826 41 


0.52004 23 


0.52184 14 


0.52439 15 


0.52714 12 


0.53584 10 


0.545907 62 



0.8051 


0.7903 


0.7779 


0.7667 


0.7337 


0.7059 24 


0.6990 


0.6446 


]l!£mcmiEEK^] 

liiinaraMeai 


0.567819 47 0.35047 


0.6061 


0.5812 


0.5557 


0.5235 


0.4894 24 


0.4296 


0.7461 19 


0.75303 77 


0.75463 95 


0.75743 89 


0.75825 81 


0.76050 84 


0.76576 81 


0.77230 90 


0.77290 70 


0.77862 49 


0.78217 90 


0.78642 86 


0.78904 46 


0.79246 79 


0.79849 75 


0.80205 84 


0.8055 13 


0.8120 11 


0.81546 98 



3.14 64 


2.24 40 


2.36 55 


11 6.3 


24 7.2 


14 4.4 


17 5.5 


35 


48 


47 21 


48 21 


55 19 


78 23 


91 55 


43 18 


55 14 


30 


24 8.7 


13 4.3 


10 4.8 


8.7 3.8 


8.6 3.7 


31 28 


5.28 

250 

5.285 

300 

5.29 

250 

5.295 

250 

5.3 

250 

5.305 

8640 

5.31 

5366 

5.3125 

2758 

5.315 

8798 

5.3175 

550 

5.32 

6178 

5.325 

200 

5.33 

300 

5.34 

350 

5.35 

400 


0.48665 


0.48918 


0.49154 


0.49334 


0.49526 


0.498846 93 


0.50211 


0.50476 


0.50775 


0.50932 


0.51268 16 


0.51513 


0.51894 


0.52239 


0.524922 98 


0.48735 25 


0.48984 14 


0.49189 25 


0.49418 54 


0.49616 27 


0.49999 10 


0.50351 17 


0.50649 18 


0.50975 21 


0.51142 41 


0.51516 18 


0.51791 36 


0.52218 20 


0.52557 33 


0.52834 10 


0.8082 


0.7916 


0.7695 


0.7510 


0.7276 


0.6983 11 


0.6598 


0.6238 


0.5830 


0.5661 


0.5144 


0.4769 


0.4298 


0.3933 


0.3677 15 



0.77044 10 


0.77587 15 


0.78178 20 


0.78810 29 


0.78906 47 


0.79808 20 


36 9.4 




4.3 1.8 


5.4 2.7 


3.0 1.2 


7.9 5.2 


17 


22 


31 7.7 


48 18 


33 


76 28 


27 26 


45 15 


15 


10 8.3 


22 


7.3 4.3 


5.24 

350 

5.26 

293 

5.27 

1434 

5.2715 

2323 

5.272 

6300 

5.2725 

3790 

5.2728 

3175 

5.2731 

2605 

5.27375 

1060 

5.275 

494 


0.479565 70 


0.48804 


0.49642 


0.50038 


0.50173 


0.50301 


0.50422 


0.50572 


0.50511 34 


giiiTiy/flrfigM 


0.480086 81 


0.48901 10 


0.49812 28 


0.50252 52 


0.50405 25 


0.50548 29 


0.50685 21 


0.50854 20 


0.50792 41 


0.51059 27 


0.7687 


0.6841 


0.5523 


0.4783 


0.4539 


0.4287 


0.4019 


0.3786 


0.3859 65 


0.7631 10 


0.7719 11 


0.78850 48 


0.79708 59 


0.79952 33 


0.80268 34 


0.80561 27 


0.80852 31 


0.80829 42 


0.81228 65 



3.4 1.6 


3.4 1.7 


68 62 


190 230 


147 90 


104 


86 


103 


100 130 


62 82 


33 






















































































































































































































































































































































































TABLE IV: (continued) 



# Traj. 

5.28 

335 

5.285 

277 

5.29 

290 

5.32 

380 

5.4 

295 


Pa 


0.511175 49 


0.51334 


0.51485 


0.522467 59 


0.537482 40 


0.51680 15 


0.51841 15 




0.2679 


0.2325 


0.2097 


{'ipDoii!) 


0.81864 63 


0.82238 67 


0.82267 54 


'^int {Per 


Tint ((V’V')) 


8.1 5.0 


14 13 


13 12 


0.526117(77) 

10.14394 


0.82835(42) 

IBfBI 

9.9(6.6 

0.541084 38 

10.09389(20) 

0.82998 44 

1 2.5(1.0 

II 

3.5 1.9 


5.3 

850 

5.325 

350 

5.35 

350 

5.375 

350 

5.3875 

350 

5.4 

350 

5.4125 

350 

5.425 

850 

5.4375 

850 

5.45 

800 

5.475 

3540 

5.475 

750 

5.4875 

3530 

5.4875 

700 

5.49375 

3390 

5.5 

2440 

5.5 

600 

5.50625 

3370 

5.5125 

3460 

5.5125 

700 

5.525 

3550 

5.525 

700 

5.5375 

3430 

5.5375 

700 

5.55 

1010 

5.55 

700 

5.575 

700 

5.6 

800 


0.46898 


0.47466 



0.48852 


0.49244 


0.49571 


0.49845 


0.50247 


0.50662 29 


0.51059 


0.51922 18 


0.51906 


0.52352 17 


0.52333 51 


0.52652 


0.52830 


0.52814 


0.53098 16 


0.53436 


0.53450 41 


0.53894 18 


0.54025 


0.542920 92 


0.54369 


0.54635 


0.54622 


0.55133 13 


0.555811 97 


0.46898 13 


0.47442 23 


0.48142 23 


0.48857 31 


0.49255 43 


0.49575 33 


0.49852 42 


0.50249 24 


0.50682 29 


0.51091 13 


0.51976 21 


0.51934 29 


0.52394 17 


0.52363 57 


0.52721 17 


0.52916 21 


0.52902 28 


0.53192 19 


0.53573 28 


0.53582 54 


0.54074 19 


0.54216 24 


0.54587 26 


0.54804 22 


0.54839 23 


0.55365 16 


0.55827 12 


1.24637 


1.23297 


1.21696 78 


1.19882 


1.18851 


1.18102 


1.17334 


1.16307 


1.15084 62 


1.14039 


1.11176 58 


1.11220 


1.09638 59 


1.0998 


1.08530 


1.07792 


1.0775 


1.06766 83 


1.0530 


1.03237 


1.02766 


1.01947 


1.01300 


1.00677 


1.00761 


0.99023 37 


0.97684 



0.94 


1.11 


1.39 37 


1.67 


3.4 1.3 


1.45 


0.91 


2.87 


3.49 91 


2.09 


8.2 1.6 


5.1 1.6 


6.4 1.1 


26 19 


6.9 1.2 


7.1 1.5 


9.8 4.4 


10 2.2 


20 6.2 


24 17 


8.4 1.6 


5.6 1.8 


3.41 


5.8 1.9 


3.20 


4.2 1.2 


1.69 30 


1.73 31 


5.3 

2050 

5.36 

16640 

5.38 

13550 

5.39 

9000 

5.4 

10000 

5.405 

8550 

5.41 

6200 

5.415 

8800 

5.42 

12300 

5.425 

14800 

5.43 

9000 

5.44 

7300 

5.46 

7350 

5.47 

1200 

5.48 

1250 

5.49 

1200 

5.5 

1300 

5.52 

900 


0.479369 57 


0.496833 32 


0.503516 45 


0.507338 76 


0.511931 89 


0.514155 98 


0.51756 


0.52114 12 


0.52368 


0.526255 80 


0.528538 79 


0.532536 72 


0.537987 39 


0.540548 88 


0.542541 75 


0.544427 63 


0.546411 75 


0.549775 69 


0.479526 69 


0.497128 34 


0.503963 47 


0.507820 82 


0.51270 10 


0.51508 10 


0.51875 20 


0.52265 14 


0.52543 13 


0.528209 93 


0.530666 88 


0.534865 86 


0.540635 40 


0.543324 97 


0.545410 79 


0.547232 66 


0.549302 69 


0.552711 78 


1.0705 


0.99318 


0.96041 


0.93911 


0.91145 


0.8973 


0.8711 


0.8458 11 


0.82691 


0.80812 


0.79081 


0.76518 


0.73047 


0.7165 


0.7047 


0.7156 10 


0.72174 33 


0.72496 37 


0.72802 47 


0.73379 45 


0.73680 53 


0.74304 61 


0.75012 48 


0.75409 39 


0.75826 36 


0.76213 39 


0.76726 50 


0.77395 45 


0.77842 93 


0.78019 95 



0.78193(90) 

4.14(96 

II 

II 

9.0(4.7 

1 0.6736 

[fDll 

0.78312 98) 

3.08 75 

8.6 4.5 


3.83 68 


7.95 70 


38 10 


6.83 85 


II 


4.20 99 


3.48 76 


12 4.2 


12 1.7 


15 2.4 


22 4.8 


30 7.0 


28 7.3 


64 


42 13 


49 12 


37 7.7 


25 5.4 


26 7.2 


15 3.4 


8.9 4.8 


11 5.4 


12 5.5 

















































































































































































































































































































































































TABLE IV: (continued) 


p 

# Traj. 

5.3 

400 

5.32 

350 

5.34 

750 

5.345 

750 

5.35 

250 

5.355 

700 

5.36 

8140 

5.365 

8565 

5.37 

5574 

5.375 

800 

5.38 

350 

5.385 

300 

5.4 

400 

5.42 

900 


5.267 


5.269 


5.271 


5.272 


Pa 


0.486532 88 


0.49252 


0.50022 


0.50303 


0.50594 15 


0.51208 


0.51635 12 


0.51993 12 


0.52249 


0.52401 


0.52644 


0.53095 


0.536151 62 


Pr 


0.48694 10 


0.49281 12 


0.50092 20 


0.50394 23 


0.50699 17 


0.50933 26 


0.51359 20 


0.51829 13 


0.52223 13 


0.52486 17 


0.52640 30 


0.52918 22 


0.53410 10 


0.539203 52 




0.92367 


0.89110 


0.8372 


0.8135 


0.7839 25 


0.7305 


0.6858 


0.6488 15 


0.6248 


0.6104 


0.5879 


0.5493 


0.51330 51 


{'lpDoli>) {Tint [Pa 



0.76815 25 


0.77747 19 


0.78332 15 43 15 



Tint [W'4’)) 


1.09 25 


1.73 54 


16 9.5 


23 ; 


11 9.3 


35 31 


77 


56 18 


52 20 


54 60 


20 


15 15 


8.3 


5.4 1.7 


5.26 

149 

5.28 

650 

5.285 

378 

5.29 

1200 

5.29125 

3960 

5.2925 

3450 

5.2925 

3118 

5.29375 

950 

5.295 

345 

5.3 

465 

5.305 

235 

5.31 

300 

5.32 

265 

5.34 

285 


0.48400 


0.492092 63 


0.49588 


0.50050 12 


0.50398 


0.504096 99 


0.50426 


0.50641 


0.50941 18 


0.512593 80 


0.515550 85 


0.517129 79 


0.520284 41 


0.525526 40 


0.48464 20 


0.493055 77 


0.49716 14 


0.50216 15 


0.50605 21 


0.50616 11 


0.50635 41 


0.50880 15 


0.515578 86 


__ ■ 


0.7814 


0.7023 10 


0.6539 


0.5884 18 


0.5299 


0.5301 16 


0.5273 


0.4896 


0.4374 41 


0.3892 


0.52044 10 



HI 

0.28789 

HI 

0.24220 


0.7605 10 


0.76843 57 


0.77561 65 


0.78489 37 


0.79314 34 


0.79319 24 


0.79334 42 


0.79804 42 


0.80598 68 


0.81064 48 


0.81585 54 


0.81816 62 


0.82037 46 


0.82483 47 


9600 


6350 


6500 


12400 


0.49225 


0.49564 


0.49997 


0.50248 63 


0.49348 28 


0.49725 54 


0.50214 54 


0.50489 72 


) 

0.6246 

) 

0.5688 

) 

0.4872 


5.273 

6800 

0.50188 

[1^11 

0.50421(79 

) 

0.449 

5.274 

9350 

0.50519 

[1^11 

0.50795(58' 

) 

0.388 

5.276 

2700 

0.50718 

mm 

0.51010(53' 

) 

0.343 


0.50995 


0.51161 


0.51318 19 


0.51496 15 


5.3 

1310 

5.31 

2200 

5.312 

1990 

5.314 

2620 

5.316 

2520 

5.318 

2600 

5.32 

2260 

5.325 

1470 

5.33 

1420 

5.34 

1000 

5.35 

900 


0.49537 


0.50213 


0.50368 


0.50495 17 


0.50750 


0.51024 


0.51128 


0.51599 


0.51875 


0.522071 90 


0.525419 69 


0.49619 16 


0.50355 23 


0.50518 27 


0.50663 18 


0.50945 44 


0.51249 25 


0.51360 31 


0.51879 18 


0.52173 10 


0.52527 10 


0.434 13 


0.2863 


0.2562 


0.77948 55 


0.78565 66 


0.79538 71 


0.80179 71 


0.79947 79 


0.80678 64 


0.81174 81 


0.81673 55 


0.81963 53 



0.7421 

KBII 

0.76363(20) 


0.6645 

[1^11 

0.77583(24) 


0.6465 

[1^11 

0.77825(24) 


0.6310(27) 1 

0.78066(26) 


0.5953 

[1^1 

0.78634(49) 


0.5554 

[1^1 

0.79230(30) 


0.5430 

mm 

0.79402(31) 


0.4773 

[1^1 

0.80338(18) 


0.4426 

mm 

0.80829(15) 


0.4049 

KSl 

0.81235(17) 


0.37128 93 0.81652 17 


8.1 3.7 


15 13 


34 23 


23 24 


7.3 5.3 


II 


2.36 90 


II 


7.2 5.9 


4.8 1.7 


9.8 6.2 


32 21 


98 62 


24 8.4 


210 220 


19 


45 68 


11 7.1 


6.8 4.9 


10 8.9 


8.6 5.5 


7.4 4.9 


112 49 


no 


no 57 


250 140 


130 


66 22 


152 


112 


240 140 


170 no 


200 120 


72 49 


42 17 


35 13 



20 10 


42 23 


57 


42 21 


160 170 


63 39 


92 75 


30 17 


20 9.5 


26 16 


11 4.9 









































































































































































































































































































































































TABLE V: Raw data from MC simulation of susceptibilities entering the specific heat. 



5.3 

1300 

5.325 

1300 

5.35 

1300 

5.375 

2930 

5.375 

1250 

5.3875 

2860 

5.3875 

1200 

5.4 

2920 

5.4 

1200 

5.40625 

2890 

5.40625 

600 

5.4125 

2890 

5.4125 

1200 

5.41875 

2900 

5.41875 

840 

5.425 

2970 

5.425 

1200 

5.43125 

2910 

5.43125 

540 

5.4375 

1300 

5.45 

1300 

5.475 

1300 

5.5 

1300 

5.525 

1350 

5.55 

1350 

5.575 

1350 

5.6 

1350 


5.335 


5.3375 


5.34 


5.345 


5.3475 


5.35 


5.3525 


5.355 


5.3575 


5.36 


5.365 


5.37 


0.0584 


0.0570 


0.0638 


0.0598 27 


0.0530 


0.0746 43 


0.0711 


0.0877 


0.0785 


0.1094 83 


0.067 


0.134 


0.0848 


0.1128 


0.105 13 


0.0926 55 


0.0824 


0.0811 49 


0.0614 


0.0521 


0.0579 


0.0425 


0.0391 


0.0379 19 


0.0338 


0.0320 16 


0.0300 


0.0407 40 


0.0375 35 


0.0416 38 


0.0321 31 


0.0588 44 


0.0732 58 


0.0626 95 


0.100 10 


0.050 10 


0.124 14 


0.069 10 


0.100 10 


0.090 13 


0.0812 63 


0.0654 94 


0.0662 49 


0.0483 96 




0.0619 44 


0.0559 38 


0.0579 36 


0.0625 32 


0.0510 34 


0.0817 46 


0.0962 61 


0.084 


0.129 12 


0.0705 97 


0.155 


0.095 


0.125 


0.112 15 


0.1089 73 


0.085 


0.0899 56 


0.069 11 


0.0532 34 


0.0622 61 


1.174 30 0.0424 78 0.0741 88 


II 


1.316 34 0.086 10 0.121 12 


1.508 41 


lIWKirfEllI 


0.209 25 


1.640 52 0.213 31 0.269 34 


1.477 42 0.155 21 0.207 24 


1.370 36 0.126 13 0.173 17 


1.221 32 0.0654 92 0.078 10 



0.0227 18 


0.0222 17 


0.0152 11 


0.0159 11 


0.0393 22 


0.0415 21 


0.0327 15 


0.0328 14 



5.3 

670 

10.0627(64) 


1 2.09(47 

5.33 

4840 

0.0732 

im 

0.0573(33) 

0.0808(36) 


5.3325 

3600 

0.0835 

0.0678(55) 

0.0911(60) 

1 2.36(26 



4000 


4000 


4000 


4000 


4000 


7130 


14800 


4000 


3800 


13200 


4000 


4000 


3100 


1200 


1300 


1300 


775 


0.0939 


0.0790 


0.114 10 


0.1090 


0.126 


0.1415 


0.1549 


0.152 


0.1086 


0.1161 


0.0806 62 


0.0750 


0.0552 


0.0483 


0.0395 21 


0.0321 


0.0264 


0.0660 58 


0.102 11 


0.0980 96 


0.120 13 


0.135 10 


0.150 10 


0.148 17 


0.1002 88 


0.1079 71 


0.0706 72 


0.0615 45 


0.0942 65 


0.130(11 

0.126 

m 

m 

0.152 

0.169(11 

0.186 

[fa 

m 

IQ 

0.182 

0.129 


0.0978 84 


0.0847 52 


2.26 22 


1.82 20 


1.98 21 


2.26 22 


2.40 18 


2.40 13 


1.94 21 


1.93 10 


1.74 16 


1.58 15 


HI 


0.0512 35 


0.0323 30 


0.0213 17 


0.0142 11 0.0293 14 


1.43 29 


1.19 18 


0.81 21 


1.19 30 


IB 

I 


0.027 64 


0.040 


0.051 


0.122 


0.073 


0.166 43 


0.105 


0.231 


0.283 


0.257 


0.263 


0.182 


0.195 


0.104 30 


0.083 


0.074 


0.053 


0.040 


0.003 


0.030 


IE 


0.019 84 


0.070 26 


0.069 32 


0.140 41 


0.104 24 


0.216 39 


0.137 35 


0.273 47 


0.363 43 


0.335 35 


0.321 63 


0.218 36 


0.243 27 


0.156 32 


0.123 31 


0.110 29 


0.038 32 


0.036 32 


0.039 23 


0.000 



5.28 

250 

5.285 

300 

5.29 

250 

5.295 

250 


0.069 


0.0531 76 0.0307 78 


0.071 11 0.052 10 


0.100 39 


nmrimni 


0.073 11 


0.126 40 
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TABLE V: (continued) 





5.3 

250 

5.305 

8640 

5.31 

5366 

5.3125 

2758 

5.315 

8798 

5.3175 

550 

5.32 

6178 

5.325 

200 

5.33 

300 

5.34 

350 

5.35 

400 


Xe.aa 


0.073 


0.1101 


0.151 11 


0.129 14 


0.195 


0.110 


0.177 


0.0579 


0.0603 


0.0724 


0.0444 57 


HI 

mi 

HI 

ni 

ni 

ni 

ni 

HI 

SI 


5.24 

350 

5.26 

293 

5.27 

1434 

5.2715 

2323 

5.272 

6300 

5.2725 

3790 

5.2728 

3175 

5.2731 

2605 

5.27375 

1060 

5.275 

494 

5.28 

335 

5.285 

277 

5.29 

290 

5.32 

380 

5.4 

295 


0.0639 67 


0.0551 


0.205 


0.46 11 


0.365 


0.371 


0.228 


0.184 


0.164 


0.123 27 


0.0479 55 


0.0597 


0.065 


0.0567 


0.0389 


Xe.cTT 


0.054 10 


0.0982 60 


0.146 12 


0.116 16 


0.196 19 


0.096 20 


0.175 14 


0.041 14 


0.0460 92 


0.060 11 


0.0262 49 


Vs—32 


0.0455 78 0 


0.195 43 


0.50 13 


0.394 54 


0.398 68 


0.231 32 


0.184 31 


0.168 49 


0.114 28 


0.0318 56 


0.050 10 


0.055 14 


0.0369 58 


0.0199 28 


Xe.TT 


.075 


1256 67 


.181 14 


.144 18 


.239 


.120 21 


.213 


.065 


.067 


.080 


0462 48 


3.311 50 


3.383 66 


3.056 87 


3.440 82 


2.50 20 


2.869 64 


0.144 



0.195 36 


0.390 


0.017 


0.316 



0.197 15 


0.316 30 


0.252 44 


0.481 51 


0.020 


0.403 38 



0642 93 



5.3 

850 

5.325 

350 

5.35 

350 

5.375 

350 

5.3875 

350 

5.4 

350 

5.4125 

350 

5.425 

850 

5.4375 

850 

5.45 

800 

5.475 

3540 

5.475 

750 

5.4875 

3530 

5.4875 

700 

5.49375 

3390 

5.5 

2440 

5.5 

600 

5.50625 

3370 

5.5125 

3460 

5.5125 

700 

5.525 

3550 

5.525 

700 

5.5375 

3430 

5.5375 

700 

5.55 

1010 

5.55 

700 


0.0468 


0.0577 


0.0552 


0.0483 59 


0.069 


0.0495 


0.0536 


0.0611 


0.0626 


0.0422 


0.0816 


0.0573 45 


0.0960 


0.084 


0.0889 


0.0902 


0.0652 


0.0846 


0.0934 


0.0703 


0.0693 46 


0.0463 


0.0469 


0.0412 31 


0.0451 33 


0.0449 


0.0248 25 


0.0342 59 


0.0336 50 


0.0290 60 


0.049 11 


0.0323 75 


0.0382 61 


0.0423 51 


0.0419 68 


0.0259 27 


0.0655 61 


0.0371 52 


0.0817 58 


mm 

0.072 10 

) 

0.094 


0.0747 56 


0.0482 73 


0.0729 52 


0.0833 82 


0.0546 48 


0.0338 66 


0.0300 33 


0.0324 37 


0.0317 39 


0.469 


0.471 


0.274 


0.226 


0.212 


0.145 31 


0.0517 56 


0.080 


0.080 


0.0591 71 


0.0351 29 


0.0447 30 


0.0552 57 


0.0483 52 


0.0545 64 


0.067 11 


0.058 


0.0647 60 


0.0606 73 


0.0453 29 


0.0871 66 


0.0558 61 


0.1031 62 


0.0957 56 


0.1025 78 


0.0655 76 


0.1092 87 


0.0729 50 


0.0529 71 


0.0510 21 


0.0536 43 


0.0568 51 


0.0513 45 


arriQ 

=0 

10(1.4) 


10(1.7) 


8.16(75 

r‘ 

9.87(74 

r‘ 

8.59 36 


raroERMl 

6.23(32 


6.50(37 

r‘ 

5.62(54 

r‘ 

6.64(81 


3.95(68 

II 

3.79(59 

r 

2.22(34 


1.80(29 

II 


0.03 10 


0.06 10 


0.38 11 


0.70 27 


0.88 13 


0.83 16 


0.423 


0.322 


0.34 16 


0.234 97 


0.002 53 


0.067 


0.057 


0.006 


0.049 


0.133 81 


0.157 89 


0.43 12 


0.77 28 


1.04 15 


0.498 92 


0.435 87 


0.41 16 


0.38 11 


0.119 65 


0.085 68 


0.093 77 


0.019 35 


0.015 29 
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TABLE V: (continued) 




5.575 

700 

5.6 

800 


Xe,(T(T 

1 Xe.CFT 

1 Xe.TT 

1 Xe, 


0.0342 


0.0324 



5.3 

2050 

5.36 

16640 

5.38 

13550 

5.39 

9000 

5.4 

10000 

5.405 

8550 

5.41 

6200 

5.415 

8800 

5.42 

12300 

5.425 

14800 

5.43 

9000 

5.44 

7300 

5.46 

7350 

5.47 

1200 

5.48 

1250 

5.49 

1200 

5.5 

1300 

5.52 

900 


0.0484 20 


0.0610 


0.0702 17 


0.0766 29 


0.0860 


0.0887 


0.1156 


0.1012 


0.1043 


0.0880 33 


0.0759 


0.0663 


0.0470 13 


0.0440 


0.0409 


0.0355 


0.0426 23 


0.0334 


0.0277 21 


0.0421 11 


imsmma 


HI 


5.3 

400 

5.32 

350 

5.34 

750 

5.345 

750 

5.35 

250 

5.355 

700 

5.36 

8140 

5.365 

8565 

5.37 

5574 

5.375 

800 


350 

5.385 

300 

5.4 

400 

5.42 

900 


0.0383 39 


0.0552 


0.086 


gl¥iyi!TIfM 


0.0638 


0.078 


0.166 


0.1316 


0.1054 


0.0667 83 


0.0696 


0.0490 


0.0435 


0.0400 


0.0722 38 


0.0738 38 


0.1066 97 


0.0918 66 


0.0958 60 


0.0758 37 0. 


0.0622 31 0. 


0.0523 31 


0.0313 12 


0.0285 24 0. 


0.0266 22 


0.0200 17 0. 


0.0261 20 


0.0157 17 


v ,=20 


0.0230 34 0. 


0.0299 51 


0.067 12 0 


0.063 11 0 


0.0444 84 0. 


0495 23 


0620 12 


0721 18 


0831 33 


0969 43 


.135 


1196 73 


1250 69 


1002 41 


0841 35 


0739 35 


0498 13 


0463 28 


1.47 22 


1.235 65 


1.260 73 


1.345 91 


1.410 88 


1.57 11 


1.52 13 


1.332 98 


1.289 81 


1.314 72 


1.205 81 


1.21 10 


0.972 81 


0.67 13 


0.015 31 


0.019 


immsm 

II 

II 
II 

II 
II 


0.072 


0.098 


0.174 


0.127 


0.152 


0.089 


0.065 


0391 20 


HI 


0324 22 


0.73 12 


0.66 12 


0.56 13 


am„- 




0.046 


0.001 


0.030 


0.012 


0.015 


0.086060 


0.049 23 


0.048 10 


0.065 13 


0.090 19 


0.120 17 


0.149 22 


0.215 35 


0.187 23 


0.207 19 


0.171 16 


0.120 17 


0.097 18 


0.062 14 


0.039 30 


0.027 28 


0.026 33 


0.016 24 


0.002 22 


0.069 16 


0.160 14 


0.1235 94 


0.0959 88 


0.0506 89 


0.054 10 


0.0396 86 


0.0255 40 


0.0209 22 


0477 53 


.085 


.085 13 


0632 90 


.096 


.194 


.154 10 











2.418(59) 

0.047 

1.943(50) 

0.050 

1.670(40) 

IBlIllli 


0.052 31 


0.056 21 


0.012 



5.26 

149 

5.28 

650 

5.285 

378 

5.29 

1200 

5.29125 

3960 

5.2925 

3450 

5.2925 

3118 

5.29375 

950 

5.295 

345 

5.3 

465 

5.305 

235 

5.31 

300 

5.32 

265 

5.34 

285 

u-- 

5.267 

9600 

5.269 

6350 


0.099 


0.0643 


0.081 


0.110 


0.253 40 


0.286 


0.303 


0.176 


0.102 


0.0706 


0.0482 


0.0666 90 


0.0360 37 


0.0364 33 


0.076 25 


0.0472 75 


0.066 23 


0.100 18 


0.261 44 


0.297 25 


0.31 10 


0.175 35 


0.082 22 


0.092 

1^ 

0.0706(87 

0.090 

[1^: 

[1^: 

0.133 

0.312(50) 

0.351 

ri^: 



5.6 1.2 


5.65 78 


4.71 62 


4.35 37 


5.46 29 


5.74 41 


4.56 45 


4.13 81 


0.039 


0.034 


0.219 


0.58 12 


0.589 


0.65 26 


0.35 11 


0.127 


) 0.0512(85) 

iiiiiwaiTBaiKBErKaii 

-( 

J 0.0372(81) 

0.0625(95) 

2.16(36 

■1 

■1 

■1 

■1 

-( 

J 0.0431(96) 

0.0552(89) 

3.02(51 

-( 

J 0.0204(29) 

0.0363(34) 

2.04(32 

0 

J 0.0163(27) 

0.0355(30) 

1.96(37 

-( 



0.123 64 


0.106 79 


0.286 51 


0.69 13 


0.680 79 


0.82 30 


0.41 12 


0.195 82 


0.012 45 


0.019 58 


0.046 47 


0.036 47 


0.029 38 


0.148 15 


0.242(45, 



0.259 53 


0.41(11 


0.327 64 


0.52(13 



































































































































































































































































































































































TABLE V: (continued) 



gaiBMi 

1 Xe.crcr 

1 Xe,crT 

1 Xe,TT 

Xe 


1 Xe„ 

IT 

1 Xe, 

T 1 

5.271 

6500 

1 0.310(40) 


1 0.323(43) 


1 0.380(49) 

1 

7.81( 


1 0.63(11) 

1 0.72(13) 


5.272 

12400 

0.445( 

m 


0.478( 

60) 


0.558( 


1 

8.40( 

’29) 

1 0.95( 

MU 

1.13( 

14) 


5.273 

6800 

0.377( 

m 


0.401( 

50) 


0.470( 


1 

7.74( 

’34) 

1 0.75(12) 

0.89( 

Ml 


5.274 

9350 

0.313( 

m 


0.328( 

59) 


0.386( 


1 

7.16( 

31) 

iBimniiiiii 

gl■Miln 

1 

5.276 

2700 

0.193( 

m 


0.195( 

39) 


0.235( 

EE) 

1 


m 

0.301( 

mi 

0.404(81 

1 

5.278 

4200 

0.1073 

(73 


0.0958 

w 

a 

0.1236 

(97 

|] 

4.20( 

;i9) 

0.1141 

:26) 

0.185(29 

1 

5.28 

4200 

0.0819 

('78 

a 



a 


m 

a 

3.88( 

; 2 o) 

0.089( 


0.136(28 

1 

1 



24 


amc= 


5.3 

1310 

1 0.108( 

m 


0.095( 

;i 6 ) 


0.124( 



3.13( 

12 ) 

1 0.138( 

mi 

O.I 661 

(32 

1 

5.31 

2200 


in 


0.164( 

19) 


0.198( 

21 ) 



mi 


mi 

0.3741 

(46' 

1 

5.312 

1990 

0.174( 

m 


0.171( 

26) 


0.211( 




mi 

0.306( 

mi 

0.372i 

(61 

1 

5.314 

2620 

0.156( 

m 


0.150( 

23) 


0.187( 



3.35( 

’ll) 

0.295( 

mi 

0.3791 

(63 

1 

5.316 

2520 

0.187( 

m 


0.187( 

41) 


0.230( 



3.37( 

18) 

1 0.38(12) 

1 0.48(11) 1 

5.318 

2600 

0.170( 

m 


0.165( 

21) 


0.199( 



3.10( 

10) 

0.333( 

mi 

0.411 

(55' 


5.32 

2260 

0.184( 

m 


0.183( 

32) 


0.225( 




mi 

0.335( 

mi 

0.409 

(65 


5.325 

1470 

0.128( 

;i9) 


0.117( 

20) 


0.143( 

:2i) 

1 

2.346 

(92) 

0.1891 


0.248 

(51' 


5.33 

1420 

0.0760 

[1^ 

l] 

0.0595 


11 

0.0824 

(82 

11 

1.883 

(75) 

0.043( 

mi 

0.084 

(16 


5.34 

1000 

0.0579 

(74 

l] 

0.0419 

(’78 

|] 

0.0615 

(86 

|] 

1.621 

(’71) 


mi 

0.070 

(16' 


5.35 

900 

10.0556(42 


0.0365(41 

a 

0.0561(43 

a 

1.456(68) 

mmssi 

1 0.060(12 

n 


TABLE VI: Raw data from MC simulation for chiral condensate susceptibility and termal suscep¬ 
tibility. 


p 

# Traj. 


Xt,cr 

Xt,T 

Xu 


Lt = 4: 


12 

amg=0.1535 

48 

5.3 

1300 



iBiiaigrem 


5.325 

1300 

0.806(31) 

BlliaiTfeMl 

IBlBWiMtai 


5.35 

1300 

0.863(37) 

BiiwifaiaKii 

tBSEsmimi 


5.375 

2930 

0.947(27) 

BlMtfiiatliW 

immuna 

issMmi 

5.375 

1250 

0.726(27) 

BiBiftifirai 

IBIUBWIBII 


5.3875 

2860 


BlKRMfcMl 


1-0.302(53)1 

5.3875 

1200 

1.17(11) 

BlKliRffaai 

IBKSUmi 


5.4 

2920 

1.396(67) 

BigraafciiW 

iBMiiimi 

IBiSMmil 

5.4 

1200 

1.117(91) 

BlKWaiigil 

IBlBIlMfigil 


5.40625 

2890 

1.95(17) 

BlWiWiBMl 

IBliailliiliW 

1 -1.37(18) 1 

5.40625 

600 

1.23(13) 

BlKiailWIl 

IBlBWillfiM 


5.4125 

2890 

2.31(22) 

miBiiniw 

iBimimi 

1 -1.81(23) 1 

5.4125 

1200 

1.43(20) 

BiWiflErfeUll 

lElKUUin 


5.41875 

2900 

1.87(13) 

BiMillTfaill 

IBllrfllllfctM 

1 -1.35(16) 1 

5.41875 

840 

1.51(17) 

■lliMtlH 



5.425 

2970 

1.585(97) 



1 -1.09(12) 1 

5.425 

1200 

1-34(14) 

BIlEBafiaiil 

IBIlBBIgWai 


5.43125 

2910 

1-316(75) 


iBEmmi 

1-0.697(81)1 

5.43125 

540 

0.99(17) 

BlKWlTfeMl 

IgiKBMIlll 


5.4375 

1300 

0.709(42) 

ElEMmi 

mammi 


5.45 

1300 

0.859(67) 

BlKBIlfcMl 

iBEmmi 


5.475 

1300 

0.486(30) 

BCTBWliai 

iBiMwaiMi 


5.5 

1300 

0.384(16) 

BilEMm 

IBWEHUH 


5.525 

1350 

0.368(15) 

1-0.1287(90) 

-0.1423(96) 


5.55 

1350 

0.297(12) 


-0.1037(70) 


5.575 

1350 

0.283(12) 

1-0.0885(74) 

-0.0909(69) 


5.6 

1350 

0.2442(92) 

1-0.0791(55) 

-0.0944(61) 
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TABLE VI: (continued) 


(3 # Traj. x 


mm 

c 

Xt,cr 

Xt,T 

Xt. 


5.3 

670 

1.61(34; 

5.33 

4840 

1.82(21’ 

5.3325 

3600 

2.20(25' 

5.335 

4000 

2.31(22; 

5.3375 

4000 

2.79(32' 

5.34 

4000 

3.89(37; 

5.345 

4000 

4.05(37; 

5.3475 

4000 

4.27(49' 

5.35 

7130 

5.56(48’ 

5.3525 

14800 

5.75(42' 

5.355 

4000 

5.24(58; 

5.3575 

3800 

3.74(39’ 

5.36 

13200 

4.29(32; 

5.365 

4000 

2.39(29; 

5.37 

4000 

2.26(24’ 


3100 

1.16(12’ 

5.4 

1200 

1.16(27’ 

5.45 

1300 

0.501(95 

5.5 

1300 

0.233(73 

5.6 

775 

0.120(31 


-1.09(12 


-1.24 13 


-1.11(12 

-1.14 

mm 

-1.13(15 

-1.23 

mm 

mm 

mm 

-1.48(15 

-1.67 

-1.58(12 

-1.85 

-1.53(19 

-1.67 

-0.99(12 

-1.07 

mm 

-1.162(97 

r 

-1.35(10 



-1.03 38 


-1.20 34 


-1.77 39 


-1.81 43 


-2.74 54 


-3.65 51 


-3.57 37 


-3.25 68 


-2.26 43 


-2.60 29 


-1.33 27 


-1.06 31 




5.28 

250 

5.285 

300 

5.29 

250 

5.295 

250 

5.3 

250 

5.305 

8640 

5.31 

5366 

5.3125 

2758 

5.315 

8798 

5.3175 

550 

5.32 

6178 

5.325 

200 

5.33 

300 

5.34 

350 

5.35 

400 


3.16 35 


2.92 28 


3.60 40 


5.10 91 


5.4 1.1 


5.27 24 


7.57 61 


6.63 77 


11 1.1 


4.70 52 


3.24 76 


3.02 60 


3.21 46 


1.39 17 



5.24 

350 

5.26 

293 

5.27 

1434 

5.2715 

2323 

5.272 

6300 

5.2725 

3790 

5.2728 

3175 

5.2731 

2605 

5.27375 

1060 

5.275 

494 

5.28 

335 

5.285 

277 

5.29 

290 

5.32 

380 

5.4 

295 


9.3 1.3 


9.3 1.4 


17 2.8 


50 10 


42 5.6 


44 7.0 


27 3.6 


39 10 


20 4.6 


15 4.0 


3.64 63 


5.65 93 


4.33 95 


1.15 13 


0.083 14 




-2.76 63 


-7 1.0 


-7 1.3 


-4.45 67 


-3.60 62 


-3.15 90 


-2.10 81 


-3.00 


-7 2.6 


-8 1.1 


-8 1.5 


-4.92 


-3 2.9 


-7 1.6 


-23 5.0 


-19 2.9 


-18 3.4 


-9 1.8 


-14 4.6 





-2 1.4 


-1.27 93 


0.21 38 


0.00 10 




















































































































































































































































TABLE VI: (continued) 


p 

# Traj. 

5.3 

850 

5.325 

350 

5.35 

350 

5.375 

350 

5.3875 

350 

5.4 

350 

5.4125 

350 

5.425 

850 

5.4375 

850 

5.45 

800 

5.475 

3540 

5.475 

750 

5.4875 

3530 

5.4875 

700 

5.49375 

3390 

5.5 

2440 

5.5 

600 

5.50625 

3370 

5.5125 

3460 

5.5125 

700 

5.525 

3550 

5.525 

700 

5.5375 

3430 

5.5375 

700 

5.55 

1010 

5.55 

700 

5.575 

700 

5.6 

800 


nSM 

c 

Xt,(T 

Xt,T 

Xt, 


X 


0.367 17 


0.355 24 


0.373 28 


0.359 25 


0.372 31 


0.367 26 


0.320 21 


0.366 18 


0.363 19 


0.0813 96 


E 


iwani 

USE 


0.591 21 


0.387 24 


0.674 23 


0.654 24 


0.675 31 


0.430 42 


0.682 30 


0.488 45 


0.513 17 


0.352 25 


0.399 10 


0.337 16 


0.295 16 


0.221 10 


0.2125 


0.147 34 


0.327 41 


0.362 40 


0.377 48 


0.487 50 


0.552 


0.185 29 


0.1341 85 -0.1474 86 -0.001 17 


!i] 


n 



0.0641 83 -0.0866 83 


0.0642 72 -0.0863 



5.3 

2050 

0.74(13) 

5.36 

16640 

0.929(49 


13550 

1.062(69 

5.39 

9000 

1.296(95 

5.4 

10000 


5.405 

8550 

1.48(10) 

5.41 

6200 

1.94(23) 

5.415 

8800 

1.87(14) 

5.42 

12300 

1.88(11) 

5.425 

14800 

1.520(84 

5.43 

9000 

1.180(77 

5.44 

7300 

1.05(10) 

5.46 

7350 

0.657(52 

5.47 

1200 

0.55(12) 

5.48 

1250 

0.526(99 

5.49 

1200 

0.49(10) 

5.5 

1300 

0.54(10) 

5.52 

900 

■iWaiKB 


0.28 


0.014 



MBiMifaCHlgIgaiTinMi 

-1.09(17 

-1.55(29 

MBlgiWiffiGMBllBMBffiHli 

-1.25(19 

-1.18(14 

MBlMifcWfiaMBlIflaiieEMi 

-1.13 11 



0.07 17 


5.3 

400 

1.384(93 

5.32 

350 

1.136(94 

5.34 

750 

2.04(29) 

5.345 

750 

2.22(27) 

5.35 

250 

2.31(27) 

5.355 

700 

2.58(46) 

































































































































































































TABLE VI: (continued) 


/3 

# Traj. 


5.36 

8140 

5.13(39) 

5.365 

8565 

4.16(31) 

5.37 

5574 

3.37(27) 

5.375 

800 

2.00(40) 


350 

1.77(32) 

5.385 

300 

1.86(38) 

5.4 

400 

0.810(98) 

5.42 

900 

■iiigweiHi 


Xt, 

(7 

-1.591 

'14' 

1 

-1.241 

10 

1 


Xt, 

r 

-1.76( 


■ 

-1.38( 

irai 

■ 


-1.09 




5.26 

149 

7.7(1.8] 

5.28 

650 

4.38(57 

5.285 

378 

5.9(1.0] 

5.29 

1200 

6.98(84 

5.29125 

3960 

18(T0) 

5.2925 

3450 

22(1.6) 

5.2925 

3118 

23(6.9) 

5.29375 

950 

13(2.7) 

5.295 

345 

5.9(1.2) 

5.3 

465 

3.91(55 

5.305 

235 

2.28(42 

5.31 

300 

3.13(70 

5.32 

265 

1.55(25 

5.34 

285 

1.23(18 


- 1.10 


0.6(2.4) 

-1.52(93 


-1.21(24 
-4.01(70 
-4.45 35 
-4(1.6) 
-2.54(58 
-1.13 24 


-1.44(27 
-4.43 77 


-5.05 


-5 1.8 


-2.99 


-1.26 


-io(ii) 

- 10 ( 1 . 1 ) 


-11 4.3 


.259 65 0.16 36 



5.267 9600 14(1.1) -2.08(29) 

5.269 6350 23(4.0) -3.96(89) 

5.271 6500 30(4.1) -5.35(80) 

5.272 12400 48(5.1) -8(1.0) 

5.273 6800 42(5.0) -7.19(98) 

5.274 9350 38(6.2) 1.25(63) 

5.276 2700 23(4.8) -3.84(84) 

5.278 4200 11.00(97) -1.64(18) 

5.28 4200 7.84(90) -1.15(17) 


-2.30 


-4(1.0 

-5.96(9( 


-9 1.1 


-8 1.0 


1.44 72 


-4.30 


-1.90 


-1.27 


-4.85(93) 


Bi^n 


-13(2.2) 

-20(2.5) 

-17(2.6) 

-16(3.3) 




-3.57(60) 

-2.51(58) 


5.3 

1310 

4.70(48 

5.31 

2200 

8.53(86 

5.312 

1990 

^CEi 

5.314 

2620 

8.6(1.1] 

5.316 

2520 

10(2.8) 

5.318 

2600 

9.5(1.1) 

5.32 

2260 

9.8(1.4) 

5.325 

1470 

7.1(1.2) 

5.33 

1420 

3.37(37 

5.34 

1000 

2.43(40 

5.35 

900 

1.56(14 


-1.021 

16 


-1.10 

-2.031 



-2.24 

-2.051 



-2.32 

-1.951 



-2.20 

-2.531 



-2.89 

-2.211 



-2.45 

-2.361 



-2.68 

-1.631 

30 


-1.76 


-2.03(32) 

-4.81(66) 

-4.65(75) 

-5.12(85) 


-5.70(76) 
-5.48(91) 
-3.46(76) 
-1.09 22 






















































































































































































































TABLE VII: Values of Xm”” merasured at the f3 nearest to the pseudocritical coupling. This values 
was taken as a constant through the critical region and added to to obtain Xm- 


Ls 

aniq 

/3 

# Traj. 

^conn 

X.m 

12 

0.153518 

5.4125 

5600 

1.01(5) 

16 

0.075 

5.35 

5000 

1.59(10) 

20 

0.04303 

5.315 

2500 

2.12(5) 

32 

0.01335 

5.2725 

120“ 

0.2(5)“ 

12 

0.307036 

5.5 

6200 

0.55(3) 

16 

0.15 

5.41 

5000 

1.00(5) 

20 

0.08606 

5.36 

2550 

1.3(1) 

32 

0.0267 

5.2925 

o 

00 

0.4(2)“ 

24 

0.04444 

5.316 

3150 

2.02(9) 


“This quantity was measured only on a small fraction of configurations. Due to limited statistics the result¬ 
ing value for x™”" is compatible with zero. However for these cases [Lg = 32 and arriq = 0.01335, 0.0267) 
Xm™ is only a small fraction of Xm and can be safely neglected within errors. 

TABLE VIIL Pseudocritical couplings /3c extimated from the reweightes curves for Xe.o-o-- 


Ls 


Pc 


0.153518 

5.4112(18) 



5.35175(82) 


0.04303 

5.3164(11) 

32 

0.01335 

5.27180(20) 


0.307036 

5.502(10) 



5.41153(61) 


0.08606 

5.36072(77) 

32 

0.0267 

5.29250(15) 

24 

0.04444 

5.3164(18) 

16 

0.01335 

5.27168(27) 
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APPENDIX B: ESTIMATE OF THE BACKGROUND OF Cy 


No significant dependence of the backgronnd Cq from the volnme of the system is expected 
since it is an nltraviolet qnantity, while dependence on amg and f3 is expected. In order to 
estimate C'o(/3, arug) we performed a linear £t of the tails in the (3 region far from the peak 
for each different valne of arUg. The procedure used is the following: 

• hrst estimate the width W of the peak. The width that we choose to use as a reference 
is the width at 75% of the total height of the peak; 

• then for each peak eliminate those points which he in the (3 region around the maximum 
of the curve at distances smaller than n ■ W, where n is a constant; 

• £t the remaining points with a linear function and study the dependence of this htted 
background on the parameter n. 

For the whole procedure to make sense, we require that the background should £t better 
as n is increased. We require the value to be stable at sufficiently large n. This proves 
qualitatively to be the case. Since the number of points at large n is not large, we choose 
the best value of n by minimizing the reduced of fh® fif of the background. 

No dependence of the background on the bare quark mass arrig was found so we can take 
C'o(/3, arriq) = Cq{(3). We choose then to do a global £t constraining all Cy peaks at different 
masses to share the same background. The result of this £t is shown in FigJTKl The best £t 
is obtained excluding a region of width VIW. Table HXl shows the stability of the £t as the 
parameter n is in the range 6-15. 

It should be noticed that even the j3 dependence is very weak and that Co{j3) is consistent 
with a constant in the (3 range examined within the statistical errors so that the whole 
procedure described here is in practice equivalent to taking Uo(/9) a constant. 
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TABLE IX: Cq{(I) fit parameters for different values of n. The fit function is C'o(/3) = O'! + ^2/3 


n 

91 

92 

X^/d.o.f. 

d.o.f. 

3 

0.608(60) 

-0.103(11) 

8.70 

84 

4 

0.535(44) 

-0.0908(84) 

5.17 

77 

5 

0.529(48) 

-0.0887(88) 

5.19 

66 

6 

0.452(44) 

-0.0751(80) 

3.56 

58 

7 

0.444(44) 

-0.0740(80) 

3.62 

53 

8 

0.436(44) 

-0.0728(80) 

3.61 

47 

9 

0.416(44) 

-0.0693(77) 

3.19 

42 

10 

0.420(44) 

-0.0695(75) 

3.11 

37 

11 

0.401(41) 

-0.0661(77) 

2.89 

34 

12 

0.400(43) 

-0.0663(83) 

2.81 

29 

13 

0.405(53) 

-0.0668(97) 

2.99 

26 

14 

0.411(56) 

-0.068(10) 

2.95 

24 

15 

0.404(56) 

-0.067(10) 

3.15 

21 
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